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Abstract. - Scale-invariant universal crossing probabilities are studied for critical anisotropic 
systems in two dimensions. For weakly anisotropic standard percolation in a rectangular- 
shaped system, Cardy's exact formula is generalized using a length-rescaling procedure. For 
strongly anisotropic systems in 1 4- 1 dimensions, exact results are obtained for the random 
walk with absorbing boundary conditions, which can be considered as a linearized mean-field 
approximation for directed percolation. The bond and site directed percolation problem is itself 
studied numerically via Monte Carlo simulations on the diagonal square lattice with either free 
or periodic boundary conditions. A scale-invariant critical crossing probability is still obtained, 
which is a universal function of the effective aspect ratio r^s = cr where r = L/t' , z is the 
dynamical exponent and c is a non-universal amplitude. 



Introduction. - The study of universal crossing probabilities in the standard random 
percolation problem at criticality has been the subject of much interest over the last decade 
^ (see ref. for a review). At the percolation threshold, Pc, on a rectangular lattice with 
length L|| and width Lj_, the probability to have at least one cluster crossing the system in one 
of the two directions (or in both) is a scale-invariant quantity in the continuum limit where 
the lattice spacing a — > 0. Actually, the crossing probability is a function 7r(r) of the aspect 
ratio, r = L±/L\\, and depends also on the boundary conditions. It remains universal in the 
sense that it does not depend on the lattice type or whether one is considering site, bond or 
site-bond percolation. 

Following the pioneering numerical work of Langlands et al. Cardy ||^ obtained an 
analytical expression for the crossing probability between opposite sides of a rectangle. He 
combined the correspondence between percolation and the limit g — > 1 of the g-state Potts 
model [H with boundary conformal field theory techniques to calculate the crossing probability 
between two segments on the boundary of a half-plane at criticality. A conformal mapping 
could then be used to transform the half-plane result into the rectangular geometry. Cardy's 
and another crossing probability formula have been recently proven rigorously Q. Monte 
Carlo j|] and exact results have been obtained for the probability to find n incipient 
spanning clusters connecting two disjoint segments at the boundary of a finite system. More 
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recently, crossing probabilities on same-spin Ising clusters in two dimensions have been also 
considered 

Let us notice that the existence of a scale-invariant crossing probability is not a trivial 
result. In the case of ordinary percolation, for example, it is linked to the vanishing, in the 
limit g ^ 1, of the scaling dimension x{q) of a boundary condition changing operator of the 
Potts model 

In this letter, we investigate the behaviour of crossing probabilities in anisotropic critical 
systems. We first give an exact treatment of weakly anisotropic standard percolation. We 
show how a length rescaling technique, which has been used in other contexts ||l^, allows 
to generalize Cardy's exact results in this situation. Then we turn to crossing probabilities 
in strongly anisotropic systems. To our knowledge, such problems have not been considered 
before, at least with the appropriate anisotropic scaling analysis. We start with the random 
walk with absorbing boundary conditions in the continuum limit which can be considered as 
a linearized mean-field approximation [|ll| for directed percolation [|T^-|l^ and illustrates on 
an exactly solvable example the behaviour of the crossing probability in strongly anisotropic 
systems. Finally we present the results of Monte Carlo simulations for the crossing probability 
in directed percolation, a strongly anisotropic problem which has so far defied all attempts at 
an exact solution, even in two dimensions. 

Weakly anisotropic percolation. - We consider the random bond percolation on a square 
lattice with anisotropic bond occupation probabilities. The system is rectangular and the 
bond probability is p\\ {p±_) in the direction parallel to the side with length L|| {L±). The 
correlation lengths, ^|| — ^ni^"^, C-L = diverge along the critical line where t — with 

the same exponent v for both directions but different amplitudes ^|| and The isotropy can 
be restored ||l^ by rescaling the lattice parameter anisotropically so that ^||a|| = ^±a±. The 
aspect ratio r is then changed into an effective one 

^11^11 -^-11/411 C-L 

In the second expression the effective aspect ratio involves the lengths measured in units of 
the corresponding correlation lengths. 

The critical amplitude ratio c — which is non-universal, is known for the q-state 

Potts model on the square lattice In the percolation limit, g ^ 1, one obtains: 

3m sinu Pc\\ l-Pc± 

c = tan — , . . ,„ r = — = , (2) 

2 sm(7r/3-w) 1-Pc|| Pc± 

which solves the problem for a weak anisotropy. When the anisotropy axis is not parallel 
to one edge of the rectangle, the rescaling leads to an isotropic system with the shape of a 
parallelogram. 

Random walk. - Standard conformal methods do not apply in strongly anisotropic sys- 
tems. The simplest example of such a system is the random walk in 1 + 1 dimensions, which 
may be considered as a directed walk in space-time. In the continuum limit, the probability 
density P{x, t) satisfies the diffusion equation 

dP d^P 
a7 = ^^- 
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Fig. 1 - Crossing probability in the time direction for a random walk with absorbing boundary 
conditions at i = and L and a uniform initial probability density at f = 0. r — /t is the 
scale- invariant aspect ratio of the system and D = 1. 



We assume a uniform initial probability density P{x, 0) = 1/L and absorbing boundary con- 
ditions at cc = and L. The solution is then obtained as the following eigenvalue expansion: 



P{x,t) 



(2p+l) — 



exp 



On the rectangle L x t, the crossing probability in the time direction is given by: 



7rr= / P{x,t)dx 
Jo 



^o(2p+l)27r2 



exp 



(4) 



(5) 



Thus the crossing probability shown in fig. |l] is a universal function of the scaled variable 
cL^/f, product of the aspect ratio r = L'^ /t, where z = 2 is the dynamical exponent for the 
random walk, by the non-universal inverse diffusion constant c = D~^. 

In strongly anisotropic systems, such a behaviour generally follows from anisotropic scal- 
ing jl6|. Under the anisotropic scale transformation, L' = L/b, t' = t/h^, the crossing proba- 
bility behaves as: 



■Kt{L,t) = b TTt 



L t 



(6) 



For a dimensionless crossing probability, Xt^^ — 0, one obtains the scale-invariant behaviour 
mentioned above with b = L. 

The probability density for a walker starting at t = with uniform probability on the 
segment [0, L] to die at time t is given by 



dt 



p=0 



exp 



Z2 



AD 

"L2 



■^2 



exp 



(7) 



where i?2(^t) — J2p^-oo u^^^^^'^^^ is a Jacobi theta function. 

Since the linearized mean-field approximation for directed percolation |jll|,|T^ is also gov- 
erned by the diffusion equation (^, the mean-field crossing probability for directed percolation 
is given by (0) when the aspect ratio is sufficiently small. 
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Fig. 2 - (a) Crossing probability ttj in the time direction as a function of the logarithm of the 
aspect ratio r = /t for site and bond directed percolation on the diagonal square lattice with either 
periodic or free boundary conditions in the space direction. The L sites are wet at t = 0. (b) Shifting 
the curves for the site problems by the same amount, the data collapse on a single scale-invariant 
universal curve for each type of boundary conditions. 



Directed percolation. - The critical crossing probability in the time direction, ttj'', has 
been studied via Monte Carlo simulations for bond and site directed percolation on the square 
lattice, with the time axis along the diagonal direction. Space and time coordinates take 
integer and half-integer values alternatively on successive spatial rows and i > 0. There are 
two directed bonds leaving a site at (x, t) and terminating on the sites at {x ± 1/2, t + 1/2). 

Sites may be wet or dry and a directed percolation cluster is a collection of connected wet 
sites starting from some source at t = (l2|-|lj]. In the bond problem, bonds are open with 
probability p. A site is wet at time t + 1/2 when it is connected through an open bond to a 
site which was wet at time t. In the site problem, sites are occupied with probability p and 
all the bonds are open. A site is wet at time i + 1/2 when it is occupied and connected to a 
site which was wet at time t. A scale- invariant crossing probability is obtained when all the 
sites are wet in the initial state at t = 0. 

The critical crossing probability was calculated for 10^ samples with sizes L = 2^^ to 2* and 
different values of the aspect ratio, r = i^/t, ranging from r ~ 2~^ to 2^. The percolation 
thresholds pj?™*^ = 0.644700185(5), = 0.70548522(4) and the dynamical exponent z = 
1.580745(10) were taken from ref. [T^]. Either periodic boundary conditions (PBC) or free 
boundary conditions (FBC) were used at / = 1 and L. The samples were grouped into 20 
packets in order to evaluate the errors on irf^. The raw data irf^ versus Inr are shown in 
fig. 2a for the largest sizes. The statistical errors are actually smaller than the size of the 
symbols. The scale invariance is manifest and finite size corrections to the leading behaviour 
are quite weak. There is a shift in In r between the curves for site and bond percolation which 
is the same for PBC and FBC. As shown in fig. 2b, an universal curve is obtained for each type 
of boundary condition after the curves corresponding to the site problem have been shifted 
by the same amount. 

In this strongly anisotropic problem, the correlation lengths, in the space direction and 
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Fig. 3 - Non-universal amplitude ratio Cb/cg as defined in eq. (|9|) as a function of tlie inverse of tfie 
system size L. The dotted line corresponds to a cubic (quadratic) fit of the data for periodic (free) 
boundary conditions. 

^11 in the time direction, diverge as 

a-ab-pc|-^ ch = e[iiP-Pc|-^'' . (s) 

Thus, generahzing the behaviour for the weakly anisotropic system in eq. (|^), for each type 
of boundary conditions, i = FBC or PBC, one expects the crossing probabiUty nf^ to be an 
universal function fi{rg,g) of the scale-invariant ratio 



veil iuy ^ ' 

in agreement with anisotropic scaling in eq. (^. As above for weakly anisotropic percolation 
and for the random- walk, the constant c is non-universal and takes different values, Cg and Cb, 
for the site and bond directed percolation on the diagonal square lattice. 

For a given type of boundary conditions, the crossing probability 7rf is the same for couples 
of values of the aspect ratio rg and rb such that the arguments of the universal function fi{res) 
are the same for the site and bond problems. This happens when Cgrg — Cb^b. Thus, with a 
logarithmic scale, the two curves have a relative shift 



In Tg — In Tb = In 



Cb 



(10) 



This shift was estimated for the different sizes for 6 points in the ascending part of the 
the curves. Since the points for the site and bond problems correspond to slightly different 
ordinates, a direct estimation of the shift was not possible. Instead we used a polynomial 
interpolation for the curve corresponding to the bond problem and the value of the abscissa 
associated with the common value of the ordinate was determined through the tangent method. 
The results for the ratio of non-universal amplitudes as a function of the inverse size are 
shown in fig. I for PBC and FBC. They converge from above and from below to the common 
asymptotic value 



Cb 



1.281(1) , 



(11) 
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corresponding to the shift Inrg — Inrb = 0.2476(8) which was used to obtain the data collapse 
in fig. 2b. 

Final remarks. - The same universal crossing probability is obtained for any finite ho- 
mogeneous density of wet sites in the initial state at t = although the finite-size corrections 
to the leading scale-invariant behaviour is stronger at large values of r (small t) when the 
initial density decreases. With a single wet site at I = L/2 in the initial state, i.e., for an 
asymptotically vanishing initial density, the crossing probability remains a function of rcff but 
the scale-invariance is lost, nf^ decaying as L~^/" . The probability 7rf'^(ri) to find n clusters 
crossing the system in the time direction at has been also studied. Details will be given 
elsewhere. 

This work can be extended by considering the crossing probability on same-spin clusters 
in strongly anisotropic spin systems, like the three-dimensional uniaxial ANNNI model at the 
Lifshitz point, for which presumably exact results for the form of the correlation functions, 
following from a generalization of local scale invariance ||l8| , have been obtained recently jl^ . 

* * * 

The author wishes to thank Lev Shchur for a stimulating introductory talk about span- 
ning clusters. 
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